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Abstract 

We derive the analytical expression of the ground state of the Hub- 
bard model with unconstrained hopping at half filling and for arbitrary 
lattice sites. 
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One of the most difficult problems in condensed matter physics is the 
characterization of the spectrum of strongly correlated Fermi systems. For 
f-D integrable systems this can be done by Bethe ansatz. For non integrable 
systems or for higher dimensional ones, no generalizations of this method 
exist, the only resources being numerical methods. Exact diagonalizations 
of the Hubbard hamiltonian which take advantage of the conservation of the 
total number of electrons N , of the total spin in the z direction S z as well as 
of the translational invariance of the lattice, has been performed by several 
authors 0, |, g. In the case of the Hubbard model, however, one would 
like to use as basis functions the simultaneous eigenfunctions of S 2 , S z , N 
and of the operators characterizing the lattice. Basis functions which are 
invariant under SU (2), U (1) and Sf (permutation group of / objects) can be 
easily constructed in terms of Young tableaux filled with quanta. By making 
projections on a particular subgroup of the permutation group Sf one could, 
in principle, construct proper basis functions for the Hubbard model on 
arbitrary lattices. Permutation group techniques were also used to solve the 
spin 1/2 Heisenberg chain Q as well as discrete quantum systems of bosons 
||, 0] . In the present paper we consider the case of the Hubbard model on a 
complete graph (unconstrained hopping) and show how the above functions 
can lead to an analytical characterization of the ground state of the system 
valid for clusters of arbitrary sizes. In spite of the simplification introduced 
by the permutational symmetry, the model keeps all the complexity of the 
standard Hubbard model coming from the continuous symmetries and its 
solution may be useful for understanding the local properties of strongly 
correlated Fermi systems. Moreover, the lacking of exact analytical results 
for the Hubbard model in dimensions D / 1 makes the analytical solutions 
derived for this model a good test for numerical computations. 
Let us start by introducing the Hubbard Hamiltonian on a complete graph 

/ / 
H = ~ t E C U CT +C/E"n n n- C 1 ) 
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where c] ,Cj a , (a =| or j) are standard fermionic creation and anni- 
hilation operators. The unconstrained hopping in Eq.([l|) makes the sys- 
tem invariant under the action of the permutation group Sf. In addition 
to this discrete symmetry there are the known symmetries of the Hub- 
bard model i.e. the U(l) conservation of the total number of electrons 
N = iVj + iVj_ = J2j ( n i T + n j|)> the SU{2) invariance under rotations of 

the total spin S of the system, and the particle-hole symmetry. To define 
the Hilbert space of the states we denote with |3 >,|2 >,|1 >,|0 >, respec- 
tively, the doubly occupied state, the single occupied states with spin up 
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and spin down, and the vacuum. The states of the system for a lattice of 
f sites are products of single site states \niri2...rif > with rn £ {0,1,2,3}, 
the dimension of the Hilbert space being just 4/. The conservation of the 
number operator allows to decompose the Hilbert space into a direct sum 
of eigenspaces corresponding to a fixed value of N. On the other hand, the 
invariance of the Hamiltonian under Sf allows to further block diagonalize 
this matrix with respect to the irreps of this group. Since the total spin of 
the system commutes with the permutation operator one can construct the 
irreps of SU (2) in terms of Young tableaux filled with quanta. To this end 
it is convenient to introduce the quantum number 

m = N + iV T (2) 

and consider all possible partitions (mi, 777,2, ■■■,rrif) of m into / parts com- 
patible with N,N^, with m; = 0,1,2,3 and mi > 777,2 > ••• > m f (in 
eq. (|2|) we restrict to > N± since the other cases follow from up-down 
symmetry). We remark that each partition is associated with an eigen- 
state I mi, 777 2 , ..., 777/ > of S z . This obviously leads to a family of states 
(eigenmanifolds of S z ) organized into levels. Each level of the hierarchy is 
characterized by a specific value of the number of "particles" m entering the 
states and by a specific value of S z = m — 3N/2. From the above family of 
states one constructs eigenmanifolds of S 2 , S z with a given St symmetry by 
filling the quanta rrii, characterizing the states in each level, in the boxes of 
a Young tableau according to the following rules. 

i) The quanta must be not increasing when moving from left to right in each 
row or when moving down each column of a given tableau. 

ii) The quanta referring to spin up and spin down states (m, = 1,2) cannot 
appear more than once in a row. 

iii) The quanta referring to doubly occupied states or to empty states (mi = 
3, 0) cannot appear more than once in a column. 

These rules directly follow from the permutational properties of the states 
1 3 >,|2 >,|1 >,|0 > and from the symmetry property of Young tableaux. To 
characterize the eigenmanifolds of S 2 , S z with S = S z (highest spin vectors), 
one must identify the filled tableaux corresponding to the highest vectors of 
SU (2) . This can be done by noting that the application of the rising operator 
on a state associated to a given filled tableaux produces a state with the same 
Sf symmetry but with a 1 changed into a 2 (1-2 flip) i.e. a state associated 
to a tableaux of the m + 1 level. We have therefore, that the tableaux that 
satisfy the filling rules also after a 1-2 flip, are the ones for which S > S z . 
It is worth to remark that by a 1 — 2 flip two different m tableaux can be 
in general associated to the same m + 1 tableau. In this case one must 
take linear combinations <j>± = tpi ± fa of the states corresponding to the m 
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filled tableaux, to get one state with S = S z (</>+) and one with S = S z + 1 
(</>_). By construction, the basis functions so obtained are simultaneous 
eigenfunctions of N,S 2 ,S Z which span the irreps of Sf. With respect to 
these functions the hamiltonian matrix has a block diagonal form with the 
dimension of the blocks given by the number of "highest weight" tableaux 
one can fill with M quanta. The degeneracy of the corresponding eigenvalues 
is just 2S+1 times the dimension of the corresponding Sf representation. We 
remark that the above set of rules can be easily implemented on a computer 
using symbolic programs like Mathematica ||. We have used a Mathematica 
implementation of our method to study clusters of several sizes. Prom these 
studies the following conjecture about the ground state of the system can 
be made: 

Conjecture: For t, U > and for N < 3 the ground state of the system is 
always of symmetry type {/}. For f < N and for N > 3 the symmetry of 
the ground state is of type {a±,a2, ...,a r } with r denoting the integer part 
of (N + l)/2, a x = f- (JV-2), ai = 2 for i = 2, ...,r - 1 and a r equal to 2 
or to 1 depending if N is even or odd. Furthermore, in the ground state S 
has always its minimal value (S = for N even and S = 1/2 for N odd). 

This conjecture appears natural since in the ground state one expects that 
the electrons will arrange (for N < f) in a configuration which minimizes the 
number of doubly occupied states and the tableaux listed above are just the 
ones which allow this. Using this conjecture we can easily derive analytical 
results for the ground state of the system. In particular, at half-filling and 
/ even, the ground state is associated to an S = (m = §/) tableau of 
type {2, 2,. ..,2}. From the filling rules it follows that there are only four 
possible ways of filling a tableau of this type with |/ quanta i.e. the ones 
reported in Fig.l. One easily check that only three of the states associated 
to these tableaux are S = states. By denoting with fa, fa, tp3, ifii the 
states associated with the tableaux from left to right of Fig. 1, we have 
that the states fa , fa and fa + -04 & re S = singlets (note that they are 
automatically orthogonal) while fa — fa is a S = 1 triplet. From the S = 
states one can construct the 3x3 block characterizing the ground state of 
the system at half filling for arbitrary f as 






Wf(f/2 + l) (3) 
/ 

(the under diagonal elements were omitted since the matrix is symmetric), 
from which the ground state energy is readily obtained as 



i 



w _i£i±£>_(fy (4) 
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with 



Z = 9 ft 2 U - y/9 pt 2 U 2 (9t 2 - pt 2 -U 2 )-3 (U 6 + pt 6 ). (5) 

We see that in both limits U — > and / — ► oo, Eq.([|) reproduces the 
known expression Eq = —t f of the ground state energy ||. Furthermore, 
the degeneracy of Eq is just the dimension of the {2, 2, 2} representation 
i.e. 

(/)' m 

It is worth to note that although the lattice is not bipartite, the ground 
state at half filling is a singlet just as for bipartite lattices. This seems to 
indicate that this property is true for arbitrary lattices. We also remark that 
the above approach can be used to compute the ground state for arbitrary 
fillings. Thus, for example, for N = 2 the ground state energy is found by 
diagonalizing the 2*2 matrix corresponding to the symmetric states with 
S=0, this giving 

E = 2 1 - f t + U/2 - {f 2 t 2 - 2tU + ftU + (U/2) 2 ) 1/2 , (7) 
while for M = 3 is obtained by diagonalizing the 3*3 block given by 



■U-(f/2-3)t -t/2y/f(f - 2) 1/2^2/ 



U - f/2 -3ty/fj2=T J . (8) 
* * -2(f - 3)t 

It is remarkable that for arbitrary fillings N < / the ground state energy 
will always be obtained by diagonalizing a matrix of rank less or equal to 
four. 

To check our result we have compared Eq.@ with numerical diagonaliza- 
tions of H up to clusters of 16 sites (this is the limit of computability on a 
supercomputer for this model). Thus, for example, for / = 16, t = 1, U = 4, 
one obtains with Lanczos method that the ground state energy is Eq = 
-12.719 while from Eq. (|) we get Eq = -12.7229. The difference 
between these two values represents the error involved in the numerical 
computations at large values of f. 
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Figure Captions Fig.l Young tableaux of type {2, 2, ...} filled with M = |/ 
quanta at N = f. The rectangular box represents | — 2 identical boxes. 
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